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ABSTRACT  (UNCLASSIFIED) 

A 

A  study  was  performed  to  investigate  the  effect  of  a  plane  interface 
between  two  media  with  different  electromagnetic  properties,  on  the 
electromagnetic  fields.  Such  an  interface  exists  between  e.g.  an 
alr/earth  or  air/water  boundary.  This  report  presents  the  results  of 
this  study. 

The  incident  field  is  assumed  to  be  a  uniform  plane  wave.  Therefore  this 
study  adresses  the  one -dimensional  scattering  problem.  The  analysis  is 
performed  in  the  s- domain,  i.e.  using  the  Laplace  transform. 

Numerical  results  are  presented  for  the  scattering  of  a  NEMP  by  the 
earth's  surface  and  the  sea. _ —  -  - - 
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SAMENVATTING  (ONGERUBRICEERD) 

Een  studie  is  verricht  naar  de  invloed  van  een  grensvlak  tussen  twee 
media  met  verschillende  elektromagnetische  eigenschappen,  op  de 
elektromagnetische  velden.  Zo'n  grensvlak  bestaat  bijvoorbeeld  tussen 
lucht  en  aarde,  of  lucht  en  water.  Dit  rapport  geeft  de  resultaten  van 
deze  studie  weer. 

Voor  deze  studie  wordt  er  vanuit  gegaan  dat  het  invallende  veld  een 
unlforme  vlakke  golf  is.  Dus.  dit  rapport  behandelt  het  den-dimensionale 
verstrooiingsprobleem.  De  analyse  is  uitgevoerd  in  het  s-domein  (Laplace 
transformatle) . 

Numerieke  resultaten  worden  gepresenteerd  van  de  verstrooiing  van  een 
NEMP  aan  het  aard-  en  het  zee-oppervlak. 
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1  INTRODUCTION 

In  most  electromagnetic  scattering  configurations  the  influence  of  the 
earth  has  to  be  taken  into  account,  because  an  object  in  the  vicinity  of 
the  earth's  surface  is  not  only  irradiated  by  the  incident  field,  but 
also  by  the  field  reflected  by  the  earth.  Such  scattering  configurations 
occur  in  Nuclear  ElectroMagnetic  Pulse  (NEMP)  interaction  studies,  which 
are  conducted  by  the  EMP-group  of  the  Physics  and  Electronics  Laboratory 
of  TNO  (FEL-TNO) . 

Therefore  a  study  was  performed  to  investigate  the  effect  of  a  plane 
interface  between  two  media  with  different  electromagnetic  properties, 
on  the  electromagnetic  fields.  Such  an  interface  exists  between  e.g.  an 
air/earth  or  air/water  boundary.  This  report  presents  the  results  of 
this  study. 

In  Chapter  2  the  formulation  of  the  problem  is  presented  in  the 
s-domain,  i.e.  using  the  Laplace  transform,  for  horizontally  and 
vertically  polarized  incident  electromagnetic  fields. 

In  Chapter  3  some  properties  of  the  reflected  and  refracted  fields  are 
discussed,  while  some  useful  approximations  are  presented  in  Chapter  A. 
In  Chapter  5  some  numerical  results  are  presented  for  the  scattering  of 
a  Nuclear  ElectroMagnetic  Pulse  (NEMP)  by  the  earth's  surface  and  the 
sea . 

Finally,  some  conclusions  are  drawn  in  Chapter  6. 
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2  DESCRIPTION  AND  SOLUTION  OF  THE  ELECTROMAGNETIC  SCATTERING 

PROBLEM  IN  THE  S -DOMAIN 

In  section  2.1  the  scattering  configuration  is  defined,  while  in 
sections  2.2  and  2.3  the  scattering  problem  is  solved  for  vertically  and 
horizontally  polarized  fields,  respectively. 

2.1  Description  of  the  scattering  configuration 

The  one -dimensional  scattering  configuration  depicted  in  fig.  2.1  is 
considered.  To  specify  the  position  in  this  unbounded  space,  we  employ 
the  coordinates  (x,  y,  z)  with  respect  to  a  given  orthogonal 
right-handed  Cartesian  reference  frame  with  origin  0.  The  unit  vectors 
along  the  axes  are  denoted  by  and  .  Then  the  position  vector  in 

this  coordinate  system  is  given  by 

(2.1)  r  -  xi^  yiyf  ziz. 

The  time  coordinate  is  denoted  by  t,  and  the  s-domain  quantities  are 
obtained  by  applying  the  Laplace  transform  with  respect  to  time.  The 
Laplace  transform  is  defined  by 

A  ® 

(2.2)  2<f(t))  -  f(s)  -  j  e'sc  f(t)  dt, 

J  o 

and  the  inverse  Laplace  transform  is  given  by 
1  r7+i"  st  A 

(2.3)  f(t)  -  2H  I  esc  f(s)  ds, 

J7-i«0 

where  i  is  the  imaginary  unit.  The  integration  for  the  inverse  Laplace 
transform  is  performed  along  the  Bromwitch  contour,  LePage  [1],  with  7  a 
positive  real  number.  In  eq.(2.3),  7  is  chosen  such  that 
exp(-7t)  | f ( t ) |  <  M,  where  M  is  a  finite  positive  real  number. 
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Furthermore,  let  the  half- space  Z)^  be  given  by 

(2.4)  Ij  -  |  x,  y,  z  6  R  |  z  <  0  ), 

and  let  the  half- space  be  given  by 

(2.5)  Z>2  -  (  X,  y,  z  e  IR  |  z  »  0  ) . 

This  configuration  is  irradiated  by  an  incident  uniform  electromagnetic 
plane  wave.  The  source,  which  generates  the  incident  field,  is  located 
in  .  The  constitutive  constants  of  the  half-space  Z)^  are  assumed  to  be 
those  of  vacuum,  i.e.  and  fi The  half-space  2)^  is  assumed  to  consist 
of  homogeneous,  isotropic,  time -invariant,  linear  and  instantaneous  and 
locally  reacting  material.  The  material  in  is  characterized  in  its 
electromagnetic  behavior  by  the  permittivity  «,  the  permeability  n,  and 
the  conductivity  o.  Furthermore  there  are  no  sources  present  in  Z)^ . 

The  incident  uniform  electromagnetic  plane  wave  propagates  in  the 
direction  of  the  propagation  vector  The  propagation  vector  x^  is  of 
unit  length  and  makes  an  angle  <p ^  with  the  normal  n  on  the  plane 
interface  between  the  two  media,  directing  inwards  V^,  and  with 
<Pi  <  90°. 

This  configuration  allows  the  waves  to  propagate  in  two  independent 
ways,  in  which  any  incident  wave  can  be  decomposed.  To  the  horizontal 
component  of  the  incident  magnetic  field  one  can  apply  a  coordinate 
transformation  such  that  this  component  runs  parallel  to  one  of  the 
horizontal  axis.  This  component  is  referred  to  as  the  TM-mode.  The 
horizontal  component  of  the  Incident  electric  field  can  be  treated  in 
the  same  way.  The  wave  propagation  in  this  case  is  referred  to  as  the 
TE-mode.  It  can  be  proven  easily  that  any  incident  plane  wave  can  be 
written  as  a  linear  combination  of  these  two  kinds  of  wave  propagation. 
Therefore,  we  shall  solve  these  two  independent  kinds  of  wave 
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propagation  in 


the  next  sections. 


Fig.  2.1.  Electromagnetic  scattering  configuration. 

a)  vertically  polarized  field, 

b)  horizontally  polarized  field. 


2.2  Vertically  polarized  incident  field 


Let  the  incident  magnetic  field  be  a  uniform  plane  wave,  and  let  it  have 
an  y-component  only  (TM-mode).  Furthermore,  the  incident  field  hits  the 
interface  between  23^  and  at  t  -  0.  Then  the  incident  magnetic  field 
is  given  by 


(2.6) 


Hy(r,s)  -  H„(s)  e 


I  €  D 


1 


with  the  propagation  vector  given  by 


(2.7) 


Xj-  sin  'P1ix+  cos 
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In  eq.(2.6),  the  propagation  constant  7q  Is  given  by 


7  -  s  ( e  n  ) 

0  0  0 


where  denotes  the  speed  of  light  in  vacuum.  Ho(s)  denotes  the  Laplace 
transformed  wave  form  of  the  incident  magnetic  field. 

Substitution  of  in  eq.(2.6)  yields 

a.  *  -7  (x  sin  <p.+  2  cos  <p.) 

(2.9)  H1  -  H  e  °  1  1  ,  re®, 

y  o  -  1 

which  represents  a  wave  propagating  in  the  positive  x-  and  2-direction. 
The  incident  electric  field  is  given  by 


(2.10) 


E^(E.s)  “  cos  <p.  Z  H1, 
x  l  o  y 

Ai  Ai 

E  (£,s)  -  -sin  2  H  . 
z  i  o  y 


where  Z  denotes  the  wave  impedance  of  vacuum  and  is  given  by 


a  a  r  l*0  I1'2 
(2.11)  Zo-Eo/Ho-  . 


The  reflected  fields  in  ®^  propagate  in  the  direction  of  the  propagation 
vector  xr“  sin  cos  V^LZ-  Since  the  scattering  process  is 

independent  on  the  spatial  variables  (x,  y),  the  reflected  fields  in 
are  given  by 

„r.  ,  :»  *  -V*r-£)  :« '  'Vx  sin  vi  z  cos  V 

H  <£,s)  -  RHe  -  R  H  e 

y  o  o 

A  A 

E^(l.s)  -  -cos  <pt  Z^,  £  €  ®: 


E^Cl.s)  “  -sin  v ^  Z^H*, 


(2.12) 
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which  are  uniform  plane  waves.  In  eq.(2.12)  we  have  used  Che  fact  that 

the  angle  of  incidence  equals  the  angle  of  reflection.  RH  denotes  the 

TM-mode  reflection  coefficient  which  must  be  determined  from  the 

boundary  conditions  at  the  plane  interface.  RB  is  given  by  the  ratio  of 

the  reflected  and  incident  magnetic  field  at  the  plane  interface,  i.e. 

RH  -  lira  Hr  /  H1.  The  expressions  for  the  electric  fields  in  eq.(2.12) 
z^O  ^  ^ 

are  obtained  after  substitution  of  the  reflected  magnetic  field  in 
Maxwe 1 1 ' s  equations . 

From  the  one -dimensional  wave  equation  the  propagation  constant  7  in  2)^ 
is  found  to  be 

(2.13)  7  -  H7II  -  (s/i(s£  +  a))U2  , 

where 


(2.14)  1  "  V*  +  7zV 

The  transmitted  waves  in  are  in  general  non-uniform  plane  waves. 
The  fields  in  and  can  only  couple  when  7^  in  2)^  equals  7^  in  2)^ 

V  x  €  IR.  Hence,  we  must  also  have  in  2)^ 

(2.15)  7  »  7  sin  (p.  . 

x  0  1 


The  transmitted  waves  can  be  represented  as 


where  TB  denotes  the  TM-mode  transmission  coefficient,  which  must  also 
be  determined  from  the  boundary  conditions  at  the  plane  interface.  TB 
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gives  the  ratio  between  the  transmitted  and  incident  magnetic  field  at 
the  plane  interface,  l.e.  T1 


lim  H*"  /  H1. 

,  n  y  y 


Z  and  Z  denote  the  normal 
n  t 


and  transverse  wave  impedances,  respectively .  The  latter  is  also  known 
as  the  surface  impedance,  Wait  [2].  Both  are  found  from  a  direct 
substitution  of  the  transmitted  magnetic  field  of  eq.^2.16)  in  Maxwell's 
equations.  After  doing  so,  we  obtain 


s«  +  o 


(2.17) 


S£  +  a 


The  propagation  constant  in  the  z-direction  is  found  from  eqs.(2.I3)  and 
(2.14),  and  is  given  by 


(2.18)  7z  -  vy-7*,  Re(72)  a  0, 

where  7  is  given  by  eq.(2.13),  and  7x  by  eq.(2.15).  Since  the 
transmitted  wave  in  2)^  must  obey  the  radiation  condition  at  z  -)  »,  we 
must  select  that  branch  of  7^  where  the  real  part  of  7^  is  positive. 

As  mentioned  before  the  transmission  and  reflection  coefficients  follow 
from  the  boundary  conditions  for  the  tangential  electric  and  magnetic 
fields,  which  state 


lim 

H1 

+  Hr  - 

lim 

HC, 

zT'O 

y 

y 

z^O 

y 

lim 

E1 

+  Er  - 

lim 

EC. 

z^O 

X 

X 

Z4.0 

X 

(2.19) 
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After  substitution  of  the  relevant  equations  in  eq.(2.19),  we  obtain 
-  Z'-  1 


(2.20) 


Z'+  1 
2Z ' 
Z'+  1 


where 


(2.21) 


Z  cos  <p. 
0  1 


n  cos  <p 


/~2  7z 

Vri  -  sm  tp. 


and  where  f)  denotes  the  index  of  refraction  given  by 


(2.22)  <7  -  ~  *  /r^r(—  +  D  • 

Observe  that  Z^cos  ip ^  corresponds  to  the  surface  impedance  of  the  medium 
in  ,  so  that  Z’  is  the  quotient  between  the  surface  impedances  of  the 
media  in  Z>^  and  Z^  ■ 

2.3  Horizontally  polarized  incident  field 


Let  the  incident  electric  field  be  a  uniform  plane  wave,  and  let  it  have 
an  y-component  only  (TE-mode).  Again,  the  incident  field  hits  the 
interface  between  Z^  and  at  t  -  0.  Then  the  incident  electric  field 
is  represented  by 


(2.23) 


Ey(l.S) 


Eo(s)  e 


V*i^ 


r  e  Z 


1 


with  given  by  eq.(2.7),  and  given  by  eq.(2.8).  Eo(s)  denotes  the 
Laplace  transformed  wave  form  of  the  incident  electric  field.  From 
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of  the  transmitted  and  incident  electric  field,  i.e.  TE  -  lim  /  E  . 

24,0  y  y 

Y  and  Y  denote  the  normal  and  transverse  (or  surface)  wave 
n  t 

admittances,  respectively.  Both  are  found  from  a  direct  substitution  of 
the  transmitted  electric  field  of  eq.(2.27)  in  Maxwell's  equations. 
After  doing  so,  we  obtain  for  the  admittances 


(2.28) 


Yt  " 

Sfl 

7x 

Y  -  — . 
n 

SM 


Again  applying  the  boundary  conditions,  we  obtain  for  the  TE-mode 
reflection  and  transmission  coefficients 


(2.29) 


Y'-  1 
Y'+  l’ 
2Y‘ 

Y'+  l’ 


where 


(2.30) 


Y  cos  ip.  u  cos  tp. 

0  1  r  1 


Y^  /z  ..  2 

t  vq  -  sin  ip. 


Notice  that  Y'  denotes  the  ratio  between  the  surface  admittances  of  the 


media  in  D ^  and 
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3  PROPERTIES  OF  THE  FIELDS 


Some  properties  of  the  reflected  and  transmitted  waves,  which  will 
clarify  some  of  the  properties  of  wave  propagation  in  the  frequency 
domain,  are  discussed  in  this  chapter.  The  expressions  in  the  s-domain 
are  transformed  to  the  frequency  domain  because  of  the  physical 
importance  of  the  frequency  domain.  This  is  simply  done  by  replacing  s 
by  iu. 

3.1  Propagation  of  the  transmitted  fields 


First  of  all  let  us  define  the  quantities  v  and  <p^.  If  the  medium  in 
is  lossless,  the  transmitted  waves  propagate  with  speed  v  given  by 


1 

(3.1)  v  - - . 


Also  for  the  lossless  case,  the  propagation  vector  of  the  transmitted 

wave  x  -  sin  <p  i  +  cos  ip  i  makes  an  angle  ip  with  the  z-axis,  where 
-t  t~x  t~z  6  t 


(3.2) 


*t  ”  l^- 


—  sin  tp. 

2  1 


After  rewriting  the  propagation  factor  in  V ^  in  the  frequency  domain, 
and  using  eqs . (3 . 1) - (3 . 2) ,  the  propagation  factor  is  written  as 


(3.3) 


- (xy  +  zy  )  -i(xk  +  zk  ) 
x  z  x  z 


where 


r  6  V 


2 


(3.4)  k  -  -iy  -  -iy  sin  <p,  -  k  sin  i p.  . 

x  x  o  yi  o  ’ 


i 

i 


I 
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and 


(3.5) 


-17, 


J  2 

Vrt  - 


sin  <p.  k 
1  0 


-  cos 

V 


^t/1  * i- 


U)£  COS  tp 


In  eqs . (3 .4) - (3 . 5)  denotes  the  wave  number  of  vacuum  and  is  given  by 

u/Cg.  Furthermore,  the  square  root  in  eq.(3.5)  must  be  chosen  such,  that 

Im(k  )  <  0. 
z 

As  mentioned  earlier,  the  transmitted  waves  in  T)^,  for  both  TM-  and 
TE-mode,  are  non-uniform  waves  due  to  losses  of  the  medium  in  .  While 
the  transmitted  fields  propagate  in  the  positive  z-direction,  they  are 
attenuated.  This  can  be  seen  easily  when  the  wave  number  k in  V ^  is 
decomposed  in  its  real  and  imaginary  part.  After  doing  so  the 
propagation  factor  is  rewritten  as 

-i(xk  +  zk  )  zk' '  -i(xk  +zk') 
x  z  z  x  z 

(3.6)  e  -  e  e 


where  k^  denotes  the  real  and  '  the  imaginary  part  of  k^ .  Both  are 
given  by 


Planes  of  equal  amplitude  are  given  by  z  -  constant,  while  planes  of 
equal  phase  are  given  by  the  equation  xkx+  zk^  -  constant.  See  fig.  3.1. 


Fig.  3.1.  Planes  of  constant  amplitude  ( - )  and  of  constant  phase 

(...) 

When  the  angle  between  the  normal  of  the  planes  of  equal  phase  and  the 
normal  of  the  planes  of  equal  amplitude  is  denoted  by  <p  ,  the  phase 
factor  of  the  propagation  factor  can  be  written  as 


(3.8) 


-i(xk  +  zk  )  -i(xk  sincp.  +  zk  cos  <p  ) 
x  z  o  i  l  i 

e  -  e 


with 


Notice  that  <p  is  frequency-dependent.  With  this  notation  we  can 
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formulate  Snell's  law  for  lossy  media,  which  is  then  given  by 

(3.10)  k  sin  <p  -  k  sin  <p.  , 

1  l  o  i 

where  ip^  can  be  interpreted  as  the  angle  between  the  propagation  vector 
of  the  transmitted  field  and  the  positive  z-axis  (see  fig. 3.1).  Using 
eq.(3.10)  the  phase  factor  can  be  rewritten  as 

-i(xk  +  zk')  -ik  (x  sin  ip  +  z  cos  <p  ) 
ii\  x  z  1  I  1 

(3 . 11)  e  -  e 

Notice  that  if  the  conductivity  a  vanishes,  < equals  <p and  eq.(3.10) 
results  in  Snell's  law  for  lossless  media. 

A  very  well-known  property  with  respect  to  wave  propagation  in  lossy 
media,  is  the  skin  depth.  It  denotes  that  depth  z  -  6,  for  which  the 
amplitude  has  been  attenuated  by  a  factor  e  ^  It  then  follows  from 
eq.(3.6)  that  the  skin  depth  S  is  given  by 


1 


with  k^'  given  by  eq.(3.7). 

The  attenuation  a  in  dB's  per  meter  follows  also  from  eq.(3.6)  and  is 
given  by 

(3.13)  a  -  -8.686  k;  ’ . 

3.2  Brewster  angle 

The  Brewster  angle  is  the  angle  of  incidence  for  which  there  are  no 
reflected  fields.  For  vertical  polarization,  we  observe  from  eq.(2.20) 
that  the  Brewster  angle  occurs  if  Z'  -  1,  i.e.  there  is  a  match  between 
the  surface  impedances  of  the  media  In  and  .  The  electromagnetic 
energy  is  then  completely  transported  into  Z'  is  given  by  eq.(2.21). 
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Therefore,  the  condition  to  be  satisfied  is 

2 

n  cos  <p. 

(3.14)  -  1  -  1. 

f 2  ~~2 

-  sin  !pi 

If  only  real  angles  of  incidence  are  allowed,  i.e.  the  incident  wave  is 
uniform,  this  condition  can  only  be  satisfied  when  the  imaginary  part  of 
t;2  vanishes.  This  occurs  when  the  medium  in  is  lossless.  Hence  from 
eq.(3.14)  with  a  -  0,  we  find  for  the  sine  of  the  Brewster  angle 


(3.15) 

Notice  that  the  Brewster  angle  exists  only  if 
For  the  special  case  p  -  eq.(3.15)  simplifies  to 


(3.16)  sin  <p* 

D 

H  0 

Since  >  1,  it  follows  from  eq.(3.16)  that  if  n  -  then  >  45  . 

It  can  be  shown  that  the  dielectric  reflection  coefficient  is  positive 
for  angles  of  incidence  smaller  than  the  Brewster  angle,  and  negative 
for  angles  of  incidence  larger  than  the  Brewster  angle. 

The  situation  RH  -  0  in  the  lossless  case  can  be  seen  as  a  special  case 
of  the  lossy  case,  but  then  with  the  condition  that  RH  has  a  minimum.  It 
is  then  customary  to  speak  of  the  pseudo  Brewster  angle.  The  pseudo 
Brewster  angle  is  frequency-dependent,  and  exists  only  if  the  angle  of 
incidence  is  larger  than  the  Brewster  angle. 

Fig. 3. 2  shows  the  modulus  of  the  reflection  coefficient  as  a  function  of 
the  angle  of  incidence  for  an  air/earth  and  air/sea  water  interface. 


TNO  report 


Page 

19 


Aftfto  of  tncMonco  [] 


Fig. 3. 2.  Modulus  of  the  reflection  coefficient  as  a  function  of  the 

-2 

angle  of  incidence  for  an  air/earth  c  -  30  and  a  -  10  S/m 
(clay),  with  tp^  -  79.7°,  and  air/sea  water  interface  -  80 
and  o-l  S/m,  with  <p^  _  83.6°.  u  -  2w  l00xl06  and  -  1. 

If  the  angle  of  incidence  equals  the  Brewster  angle,  the  high-frequency 
components  are  not  reflected,  because  then  the  reflection  coefficient 
vanishes.  In  those  cases  that  the  angle  of  incidence  deviates  from  the 
Brewster  angle,  the  reflection  coefficient  reduces  to  the  high-frequency 
limit.  It  will  be  shown  in  Chapter  U  that  this  limit  is  the  dielectric 
reflection  coefficient,  i.e.  the  reflection  coefficient  with  zero 
conductivity. 

Figs.  3.3a,  b  and  c,  show  some  of  the  above  mentioned  properties. 
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(C) 

Fig. 3. 3.  Reflection  coefficient  as  a  function  of  frequency.  Brewster 
angle  <p^  -  60°,  -  3,  /»  -  1,  a  -  10  2  (dry  soil). 

a)  (pi  -  50°, 

b)  <p.  -  60  , 

c)  -  70  . 

Fig. 3. 3a  shows  the  reflection  coefficient  as  a  function  of  frequency  if 
the  angle  of  incidence  is  smaller  than  the  Brewster  angle.  Fig. 3. 3b 
shows  the  reflection  coefficient  if  the  angle  of  incidence  is  equals  the 
Brewster  angle,  and  fig. 3. 3c  if  the  angle  of  incidence  is  larger  than 
the  Brewster  angle.  In  all  cases  the  Brewster  angle  -  60°. 

If  the  incident  electric  field  is  horizontal  polarized  a  similar 
phenomena  like  the  above  described  Brewster  angle  can  occur.  The 
condition  to  be  satisfied  is  then  given  by 


(3.17) 
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which  yields  for  the  Brewster  angle  for  horizontal  polarization 


(3.18) 

Notice  that  the  Brewster  angle  qp*  exists  only  if  11^  >  t  . 

In  most  practical  applications  »  1.  Consequently,  the  Brewster  angle 
for  horizontal  polarization  does  not  occur. 
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4  USEFUL  APPROXIMATIONS 


In  some  practical  cases  useful  approximations  can  be  made,  that  reduce 
the  complexity  of  the  equations  given  in  the  previous  chapters.  Two  kind 
of  approximations  can  be  made  : 

-  Low- frequency/high- loss  approximation, 

-  High- frequency/low- loss  approximation. 

These  approximations  will  be  presented  in  the  next  two  sections. 

4.1  Low-frequency/high-loss  approximation 


For  the  special  case  that  the  spectrum  of  the  incident  field  contains 
low  frequencies  only,  some  useful  approximations  can  be  made.  We  start 
with  approximating  the  index  of  refraction  given  by  eq.(2.22).  Now  if 
P  1,  i?  can  be  rewritten  as 

w  ( 

(4.d  ,^/Ifa-i). 

o 

If  — —  *>  1 ,  we  certainly  have  that  P  1  -  — - —  sin2a>.  .  Then  k  is 
we  w  e  fr^r  1  z 

approximated  by  (cf.  eq.(3.5)) 


(4.2) 


where  we  have  used  rj  given  by  eq.(4.1).  From  eq.(3.12)  together  with 
eq.(4.2),  we  find  for  the  skin  depth 


(4.3) 


S  = 


_2_ 

now 


Fig.  4.1  shows  some  typical  values  for  the  skin  depth,  with  y,  - 
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Fig  4.1.  Skin  depth  for  various  conductivities  with  p  - 
(low-frequency/high-loss  approximation) . 

From  fig.  4.1,  we  can  conclude  that  the  low-frequency  components  of  the 
transmitted  waves  in  .  can  travel  large  distances  without  significant 
attenuation. 

Also  using  the  low-frequency  approximation,  k^  is  approximated  by 

(4.4)  k  *  i)'k  , 

10 

where  q'  denotes  the  real  part  of  q  given  by  eq.(4.1).  From  Snell's  law, 
it  then  follows  that 

1 

(4.5)  sin  <pl  -  -  sin  «  0, 

q' 

which  implies  that  <p^  =  0.  From  this  we  can  conclude  that  for  low 
frequencies  the  transmitted  waves  propagate  almost  entirely  in  the 
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direction  of  the  positive  z-axis. 


The  low-frequency  approximation  can  also  be  applied  successfully  to  the 
expressions  for  the  reflection  and  transmission  coefficients.  We  do  so 
by  approximating  Z'  and  Y'  first.  After  substitution  of  eq.(4.1)  in 
eqs.(2.21)  and  (2.30),  we  obtain 


t;  cos  <f> ^ 


(4.6) 


Y'  = 


with  r\  again  given  by  eq.(4.1).  The  reflection  and  transmission 
coefficients  for  a  vertically  polarized  incident  field  then  simplify  to 


*H  =  l  C°S  VL  '  "r 

ri  COS 

(4.7) 

-  2 1)  cos 

TH  =  - — 

r)  cos  +  Mr 


For  very-low  frequencies  and  for  angles  of  incidence  not  near  the 

O 

Brewster  angle  (which  is  slightly  less  than  90  in  this  case),  Z'»>  1. 
Therefore  RH  =  1 ,  and  T8  =  2 . 

After  substitution  of  the  approximated  Y'  in  the  expressions  for  the 
reflection  and  transmission  coefficients  for  a  horizontal  polarization, 
we  obtain 


*  u  cos  <p  -  v 

RE  =  — - - - 

fircos  <Pi  +  V 


T 


E 


2  cos  <pi 
(*rcos  <Pi  +  ri 


(4.8) 
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For  low  frequencies,  ij  ►  1.  Hence,  RE  =  -1,  and  TE  =  0. 

Notice  that  with  using  the  low-frequency  approximation,  all  quantities 
have  become  independent  on  < . 

A. 2  High-frequency/low-loss  approximation 


In  the  same  way  as  in  section  4.1,  some  useful  approximations  can  be 

made  for  high  frequencies.  The  high-frequency  approximation  is  valid  if 

— —  <t  1  -  — - —  sin2<p.  .  Using  a  first  order  approximation  for  the  index 
r^r 

of  refraction,  r]  simplifies  to 


(4.9) 


Up.  (l  + 

r  r  2iwe 


The  propagation  constant  kz  is  given  by 


cos  <o 


(4.10) 


(u> 


2ie  cos  <p 


which  yields  for  the  skin  depth 


2«  c 


(4.11) 


00  / 


COS  <p 


V 


t  u  -  sin  ip.  -  2 
r  r  ’i 


Notice  that  the  skin  depth  has  become  frequency-independent. 

After  substitution  of  eq.(4.9)  in  (3.6),  we  obtain  for  the  propagation 
factor 


-i(xk  +  zk  )  -z/6  -i  -  (x  sin  ip  +  z  cos  <p  ) 

(4.12)  eXZ=eeV  E  C 


Also  using  the  high-frequency  approximation,  the  TM-mode  reflection  and 
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transmission  coefficients  can  be  written  as 


V.  -  vC 


€  cos  <p.  -  Ve  u  -  sin  c p. 
r  i  r  r  1 


(4.13) 


€  COS  <p.  +  sr, 


- _ _ _  -u-sintp. 

r  Ti  r  r  i 


2e  cos  a>. 
r  1 


ip.  +  /■ 


e  cos  ip.  +  Ve  u  -  sin  p. 
r  l  r  r  i 


l-i/ 

l+i/’ 

2 

1+t/’ 


with  i/  given  by 


(A. 14) 


i/ 


(Mr/«r) 


1/2 


COS 


cos 


The  TE-mode  reflection  and  transmission  coefficients  are  in  this  case 
given  by 


(4.15) 


LL  COS  tp .  - 

r  i 


/i  cos  tp.  + 
r  i 


J  •  2 

v (  p  -  sin  tp. 
r  r  i 


vC 


c  M  -  sin  <p, 
r  r  i 


1-p 

1+p 


2p  cos  <p. 


u  cos  tp,  + 
r  i 


/ : 


-  sinZip, 


2 

1+p’ 


with  p  given  by 

cos  (p 

(4.16)  p  -  (er//ir)1/2  - 

cos  tp^ 

Thus,  when  the  high-frequency  approximation  can  be  used,  we  can  make  the 
following  observations  for  both  time  and  frequency  domain  as  well  as  TM- 
and  TE-mode.  The  reflected  waves  correspond  to  the  lossless  case,  i.e. 
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as  if  a  -  0.  The  transmitted  waves  are  also  equal  to  those  of  the 
lossless  case,  except  that  there  is  a  position-dependent  attenuation 
with  a  factor  e  . 
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5  SCATTERING  OF  A  NEMP 

In  this  chapter  results  are  presented  for  an  incident  Nuclear 
ElectroMagnetic  Pulse  (NEMP) .  The  NEMP  is  a  transient  signal  with  a 
very-short  rise  time  of  about  5  nsec  and  a  large  electric  field  strength 
with  a  peak  value  of  50  kV/m.  Because  of  the  short  rise  time  the 
spectrum  of  the  NEMP  stretches  up  to  200  MHz. 

5.1  Definition  of  the  NEMP 


Usually,  the  NEMP  is  approximated  by  a  double  exponential  function  given 
by 

(5.1)  Eo(t)  -  A  (e'at-  e'^), 
with 

A  -  5.278  x  10*  [-] , 

m 

(5.2)  o  -  3.705  x  106  [s'1]  , 

0  -  3.908  x  108  [s'1]  . 


The  corresponding  spectrum  is  given  by 

(5.3)  E  (iu)  -  A(-~-  -  1-^3). 

o  lw+a  11V+/9 


Fig.  5.1  shows  the  wave  form  and  the  spectrum  of  the  NEMP. 
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5.2 


Scattering  of  the  NEMP  by  the  earth's  surface 


The  constitutive  constants  e  and  o  of  the  earth  depend  strongly  on  the 
structure  and  the  kind  of  the  soil.  Table  5.1  gives  an  overview  of  some 
typical  soils. 


soil 

€r 

N 

o 

[S/m] 

f  »  a/2ne 
[MHz] 

Humid  soil,  Clay 

30 

0.01-0.02 

6 

Fertile  cultivated 

15 

0.005 

6 

Grass ,  Meadow 

3.6 

0.05-0.11 

252 

Dry  soil 

4 

0.01 

45 

Table  5.1.  Constitutive  parameters  of  some  typical  soils. 

For  all  soils  listed  in  table  5.1  the  relative  permeability  ■*  1 .  The 
last  column  in  table  5.1  gives  the  frequency  f  for  which  f  »  o/2ne. 
Therefore,  the  low-frequency  approximation  is  valid  for  all  frequencies 
much  smaller  than  the  values  listed  in  this  column.  From  fig. 5. lb  one 
can  see  that  for  all  of  these  values  the  incident  NEMP  does  still 
contain  a  considerable  amount  of  energy,  so  that  the  low-frequency 
approximation  may  not  be  applied  for  the  frequencies  that  are  contained 
in  the  NEMP. 

Figs.  5. 2-5. 4  show  the  scattered  field  strength  in  the  frequency  and  the 
time  domain.  The  incident  field  is  assumed  to  be  a  horizontally 
polarized  NEMP. 


0.00 


Frequency  [HxJ 


(a) 


(b) 

Fig. 5. 3.  Transmitted  electric  field  strength  for  dry  soil,  -  4, 

o  -  10  2,  -  1 ,  at  z  -  0,  5,  10  m.  Horizontally  polarized 

incident  NEMP,  -0°. 

a)  frequency  domain  (1  -  Eo(0)  “  1.41x10  2  jjj  s  '), 

b)  time  domain  (1-50  kV/m) . 
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The  time-domain  results  were  obtained  by  applying  a  1024  points  inverse 
Fast  Fourier  Transform  (FFT) ,  with  a  sampling  frequency  of  1  MHz .  From 
fig. 5. 3  and  5.4  one  can  clearly  observe  that  the  high-frequency 
components  of  the  transmitted  fields  are  attenuated  with  increasing 
depth  z.  This  is  the  effect  of  the  skin  depth  6,  which  decreases  for 
high  frequencies. 

5.3  Scattering  of  a  NEMP  by  water 

The  constitutive  constants  of  sea  and  fresh  water  differ  considerably. 
See  table  5.2. 


€ 

r 

O 

f  -  o/2*(. 

0 

[S/m] 

[MHz] 

Sea  water 

80 

1-5 

11CP 

Fresh  water 

80 

0.001-0.1 

22 

Table  5.2.  Constitutive  parameters  of  water,  /»  -  1. 

Because  of  the  relative  large  conductivity  a  of  sea  water,  it  can  be 
considered  as  a  good  conductor  for  the  frequencies  of  interest.  This 
observation  is  supported  by  the  last  column  of  table  5.2,  which  again 
lists  the  frequency  f  given  by  f  -  o/2*re  .  Therefore,  if  the  incident 
NEMP  is  scattered  by  sea  water  the  low- frequency/high- loss  approximation 
can  be  used  for  all  frequencies  of  the  NEMP  spectrum. 


From  the  low-frequency/high-loss  approximation,  we  conclude  that  a 
horizontally  polarized  Incident  NEMP  is  reflected  almost  totally  by  sea 

0  E 

water,  with  a  phase  change  of  180  (recall  R  =  -1).  Now  the  electric 
field  strength  Et  -  ||£t||  just  beneath  the  surface,  is  given  by 


(1  +  i)  cos  E  , 


Z4O 


(5.4) 
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where  E1  denotes  the  field  strength  of  the  incident  field,  i.e. 

E*  “  l|Ei||  -  E^.  Since  o/of  >  1,  it  follows  that  the  electric  field 
strength  just  beneath  the  surface,  is  small  compared  with  the  incident 
electric  field  strength. 

The  magnetic  field  strength  H*"  -  ||Ht||  just  beneath  the  surface,  is  given 
by 


(5.5) 


^  +I^7>1/2  £t- 


Zj.0 


where  Y  denotes  t:.e  dielectric  impedance  of  the  medium  in  t> ^  given  by 
(5.6)  Y  -  (t/M)1/2. 


Using  the  relation  E*"  -  TE  E* ,  we  find  after  some  simple  manipulation 
(5.7)  HC  -  2  cos  <pi  H1. 


Hence,  the  magnetic  field  strength  can  be  twice  the  incident  magnetic 
field  strength. 

For  a  vertically  polarized  incident  NEMP,  the  reflection  coefficient 
equals  almost  1,  and  so  the  transmission  coefficient  equals  almost  2. 
Therefore,  the  magnetic  field  just  beneath  the  surface,  is  twice  the 
value  of  the  magnetic  field  strength  of  the  incident  field.  The  electric 
field  strength  EC  just  beneath  the  surface,  is  given  by 


(5.8) 


7 — — —  HC 
iu«+<7 


Z  H 


1  + 


iue 


zj.0 


where  Z  denotes  the  dielectric  impedance  of  the  medium  in  given  by 


Z  -  (ji/O 


(5.9) 
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t  i 

Using  the  approximation  H  »  2  H  ,  we  finally  arrive  at 

(5.10)  EC  -  2  e'1/2  ( - )1/2  E1  z,0 

'  1  +  A 

lue 

1 

Since  -  4  1,  ve  can  conclude  from  eq.(5.10)  that  the  electric 

1  +  T"  - 

field  strength  in  is  very  small  compared  with  the  incident  electric 
field  strength. 

From  these  arguments  we  can  conclude  that  in  the  case  of  scattering  by 
sea  water,  for  both  horizontal  and  vertical  polarization,  the  electric 
field  strength  in  is  relatively  small  compared  with  the  incident 
electric  field  strength,  while  the  magnetic  field  strength  in  ,  equals 
twice  the  incident  magnetic  field  strength.  However,  the  transmitted 
fields  are  attenuated  rapidly  for  increasing  depth  z  due  to  the  skin 
depth  6 . 

Figs.  5. 5-5. 7  show  the  scattered  field  strengths  for  vertical 
polarization,  both  in  the  time  and  frequency  domain.  The  time-domain 
results  were  obtained  using  a  1024  points  inverse  FFT.  The  sampling 
frequency  was  1  MHz . 
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(a) 
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Omt  [mk] 

(b) 

Fig. 5. 6.  Transmitted  magnetic  field  strength  from  sea  water  (e 
a  -  1,  n  -  1) .  at  z  -  0,  1,  5m.  Vertically  polarized 
incident  NEMP,  < p.  -  0°  . 

1  A  , 

a)  frequency  domain  (1  -  YgE^O)  -  3.74x10  5  —  s  *)  , 

b)  time  domain  (  1  -  132  A/m). 


Transmitted  electric  field  strength  from  sea  water 
a  -  1,  /»  ■*  1)  p  at  z  -  0,  1,  5  m.  Vertically  polarized 

incident  NEMP,  -  0°. 

a)  frequency  domain  (1  -  E^CO)  “  1 • 41x10  2  ^  s  l) , 

b)  time  domain  (  1  «  50  kV/m) . 
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6  CONCLUSIONS 

The  expressions  for  the  scattered  electromagnetic  fields  can  be 
approximated  using  the  low-frequency/high-loss  or  high-frequency/low- 
loss  approximation.  It  was  shown  that  the  spectrum  of  the  NEMP  is  such 
that  only  in  the  case  of  scattering  by  sea  water  the  low-frequency/ 
high-loss  approximation  can  be  used.  In  all  other  cases  no  approxima¬ 
tions  are  allowed,  that  is  if  the  exact  solution  of  the  scattering 
process  is  searched  for. 

It  was  also  shown  that  in  the  case  of  scattering  by  sea  water  and  if  the 
amplitude  spectrum  of  the  incident  field  is  such  that  the  low-frequencv/ 
high-loss  approximation  may  be  applied,  the  transmitted  magnetic  field 
strength  can  be  twice  the  value  of  the  incident  field  strength.  In  this 
case  the  transmitted  electric  field  strength  is  very  small  compared  to 
the  incident  electric  field  strength. 

The  analysis  and  results  of  the  previous  chapters  where  performed  in  the 
Laplace  or  frequency  domain.  Time -domain  results  can  be  obtained  by 
applying  the  inverse  Fast  Fourier  Transform  (FFT)  numerically  to  the 
expressions  for  the  reflected  and  transmitted  waves. 

Although  the  expressions  for  the  scattered  electromagnetic  fields  are 
rather  complicated,  it  is  believed  that  time -domain  expressions  can  be 
arrived.  This  is  a  subject  of  further  research. 


Ir  R.  Middelkoop  Ir  J.J.A.  Klaasen 

(EMP  group  leader)  (author) 
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